Self-affine fractality in 7r + p and K + p collisions at 250 GeV/c 
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Abstract: Taking into account the anisotropy of phase space in multiparticle production, a 
self-affine analysis of factorial moments was carried out on the NA22 data for 7r + p and K + p 
collisions at 250 GeV/c. Within the transverse plane, the Hurst exponents measuring the 
anisotropy are consistent with unit value (i.e. no anisotropy). They are, however, only half 
that value when the longitudinal direction is compared to the transverse ones. Fractality, 
indeed, turns out to be self-affine rather than self-similar in multiparticle production. In 
three-dimensional phase space, power-law scaling is observed to be better realized in self- 
affine than in self-similar analysis. 
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1 Introduction 



The suggestion that normalized factorial moments F q (5y) of particle- multiplicity distributions 
in ever smaller phase-space intervals Sy, may show power-law behavior [1] 

F q (5y)(K(Sy)-^ , (Sy - 0) (1) 

has spurred a vigorous experimental search for linear dependence of lni^ on — hi Sy [2]. 
Analogous to a similar behavior at the onset of turbulence, such a dependence is referred to 
as "intermittency" . Power-law scaling is typical for fractals [3], i.e., for self-similar objects 
of non-integer dimension. The powers (ft q are related to the anomalous dimensions d q = 
4> q /(q — 1) measuring the fractality of a system [4]. In general, however, only approximate 
scaling has been observed in the experimental data. 

When comparing log-log plots for one phase-space dimension, one notices that the \nF q 
saturate at small Sy. This can be explained as a projection effect of a three-dimensional 
phenomenon [5]. In three-dimensional analysis, however, the power law also does not hold 
exactly in all data. In NA22, for example, the 3D data are seen to even bend upward [6]. 

A deviation from exact scaling can be expected from an anisotropy of occupied phase- 
space. To account for such an anisotropy, Wu and Liu have suggested [7] that the scaling 
property should be different in longitudinal and transverse directions and the local multiplicity 
fluctuations are self-affine rather than self- similar. If that is the case, the anomalous scaling 
of factorial moments can be observed to be retained exactly, only under a self-affine analysis, 
where the shrinking ratio is allowed to be different in different directions. In a self-similar 
analysis, in contrast, all directions are forced to have an identical shrinking ratio. It should 
be remembered, however, that the scaling law is expected to be distorted in full phase space 
due to correlations imposed by momentum conservation. 

The experimental data sample is described in Sect. 2. The method of self-affine analysis 
is briefly summarized in Sect. 3. The results of the self-similar analysis are shown in Sect. 4. 
After reducing the influence of momentum conservation in the full region, so-called Hurst 
exponents are obtained by fitting one-dimensional factorial moments as shown in Sect. 5. 
In Sect. 6, the results from the self-affine analysis compared to those from the self-similar 
analysis are given. Conclusions are summarized in Sect. 7. 

2 The data sample 

In the CERN experiment NA22, the European Hybrid Spectrometer (EHS) was equipped with 
the Rapid Cycling Bubble Chamber (RCBC) as an active target and exposed to a 250GeV/c 
tagged, positive, meson enriched beam. In data taking, a minimum bias interaction trigger 
was used. The details of spectrometer and trigger can be found in [8,9]. 

Charged particle tracks are reconstructed from hits in the wire- and drift-chambers of 
the two-lever-arm magnetic spectrometer and from measurements in the bubble chamber. 
The momentum resolution varies from 1-2% for tracks reconstructed in RCBC, to 1-2.5% for 
tracks reconstructed in the first lever arm and is 1.5% for tracks reconstructed in the full 
spectrometer. 

Events are accepted for the analysis when the measured and reconstructed charge multi- 
plicity are the same, charge balance is satisfied, no electron is detected among the secondary 
tracks and the number of badly reconstructed (and therefore rejected) tracks is 0. The loss 



of events during measurement and reconstruction is corrected for by applying a multiplicity- 
dependent event weight normalized to the topological cross sections given in [9]. Elastic 
events are excluded. Furthermore, an event is called single-diffractive and excluded from the 
sample if the total charge multiplicity is smaller than 8 and at least one of the positive tracks 
has a Feynman variable \xp\ > 0.88. After these cuts, the inelastic, non-single-diffractive 
sample consists of 59 200 7r + p and K + p events. 

For laboratory-momenta plab < 0.7GeV/c, the range in the bubble chamber and/or the 
change of track curvature is used for proton identification. In addition, a visual ionization 
scan is used for plab < 1.2GeV/c on the full K + p and on 62% of the 7r + p sample. Positive 
particles with plab > 150GeV/c are given the identity of the beam particle. Other particles 
with momenta plab > 1-2 GeV/c are not identified in the present analysis and are treated 
as pions. 

In spite of the electron rejection mentioned above, residual Dalitz decay and 7 conversion 
near the vertex still contribute to the two-particle correlations. Their influence on our results 
has been investigated in detail in [10]. 



3 The method 

In the language of self-affine analysis in three-dimensional phase space (here denoted as 
Pa,Pb,Pc), only under the self-affine transformation 5p a — ► 5p a /X a , 5p b — > 5pb/X b , 8p c — > 
dPc/Xc with non-identical shrinking ratios A a , X b and A c , are the factorial moments expected 
to have the well-defined scaling property 

.(a) ,(b) Ac) 

F q (5p a ,5p b ,5p c ) = xi" \t q Ac F q (X a Sp a ,X b 8p b ,X c S Pc ). (2) 

The shrinking ratios A a , X b , A c are characterized by the so-called roughness or Hurst exponents 
[11] 

In A- 

Hij = ^nX ±, = a > & or a ' c or ft ' c )' ( 3 ) 

with 

Xi < Xj, < < 1, (4) 

describing the anisotropy of the system under study. For Hij = 0, Aj = 1, the scaling property 
does not exist in direction i, only in direction j. For H^ = 1, the self-affine transformation 
reduces to a self-similar one, meaning that the system is isotropic in these two directions. 
For < H^ < 1, non-trivial self-affine fractality exists in the plane, i.e, the fluctuation 
is anisotropic in that plane. 

The Hurst exponents can be deduced from the data by fitting three one-dimensional 
second-order factorial-moment saturation curves [5] 

F 2 (i) (Mi) = Oi- PiM~^ , (i = a, b, c) (5) 

where Mj = Api/5pi is the number of sub-divisions in direction i, Apt and Spi are the initial 
and final interval size in direction i, respectively, and and 7$ are three fit parameters. 

The Hurst exponents are determined from the parameter 7^ as 

1 + li 

Hij = , (i,j = a,b or a,c or b,c). (6) 

J- + 7i 



With these Hurst exponents, a self-affine analysis can be executed according to (2). If self- 
afhne fluctuations of multiplicity do exist in multiparticle production, exact scaling should 
be observed in three-dimensional phase space. 

A scaling function similar to (2) for two variables has also been suggested by J. Wosiek 
[12] as a requirement for hyperscaling from a formal analogy with statistical physics. 

4 Self-similar analysis 

The results of a self-similar analysis in 1-, 2- and 3-dimensional phase space are presented in 
Fig. 1. The initial intervals for the three phase-space variables, rapidity y, azimuthal angle 
if and transverse momentum p t , are defined as: 

-2 < y < 2 

< if < 2vr 

0.001 < Pt < 10 GeV/c. 

To avoid trivial effects due to lack of translational invariance, all variables are transformed to 
the corresponding cumulative variables by the Ochs method [13]. The experimental resolution 
in y, (p and pt has been studied in detail in [10]. The limited available statistics, rather than 
the experimental resolution, sets the limit on the smallest bin size to be used in the analysis. 

In one-dimensional projection, the partitioning M = 1, 2, ... 40 is used for all three vari- 
ables. In two-dimensional projection, the partitioning in each direction is M = 1, 2, . . . 20, so 
that for the area it is M2 = l 2 , 2 2 , . . . 20 2 . In the three-dimensional case, M = 1, 2, ... 15 is 
used for each direction, so that M3 = l 3 , 2 3 , . . . 15 3 for a three-dimensional box. 

In the 1-D analysis (first column of Fig. 1), F2 saturates at three different values when 
using y, p t and 92, respectively. In the case of y, F2 increases rapidly with increasing InM 
at small M and reaches a saturation value which is the highest of the three. The trend is 
followed for the case of p t , but at lower values of F2. When using ip, on the other hand, F2 
increases with increasing In M only above an initial decrease. 

From the 2-D analysis (second column of Fig. 1) in the (y,pt) plane, we observe an onset 
for a saturation at medium InM, followed by an upward bending at large InM. An upward 
bending is observed in the (y, ip) plane. In the {pt, ip) plane, a decrease at low InM is followed 
by an upward bending. 

In the 3-D analysis, lnF q is bending upward for all orders of q. 

5 Hurst exponents for higher-dimensional phase space 

Momentum conservation by itself causes a correlation and can, therefore, distort the scaling 
behavior expected from the dynamics of particle production [14]. The influence of momentum 
conservation on the factorial moments is expected to be different in the various variables. The 
variable p t contains only the absolute value of momentum in the transverse plane without any 
information on the direction. The influence of momentum conservation, therefore, is small 
for this case. For rapidity y, the influence of leading particles is reduced by the y-cut given 
above. Therefore, the influence of momentum conservation in y is not significant. For the 
variable <p, however, all directions in the transverse plane are included for M = 1, so that F2 



is dominated by transverse momentum compensation, which explains the decrease of i^bC^O 
with increasing InM at low InM as shown in Fig. 1. 

After reducing the influence of momentum conservation by excluding low values of M, it 
is easy to obtain the Hurst exponents from the data by means of (5) and (6). The fit results 
obtained according to (5) are shown in Fig. 2 for all phase-space variables considered. The 
parameter values are given in Table 1. Accordingly, the Hurst exponents deduced from (6) 
are: 

flypt = = 0.48 ± 0.06 ; 



2.14 
1.0139 



H y!f = —r-r— = 0.47 ± 0.06 



2.14 



^ = I^= -" ±0 - 01 - 
From these Hurst exponents, we, indeed, observe anisotropy (H y j « 0.5) between longitu- 
dinal and transverse directions of multiparticle production, while there is an isotropy in the 
transverse plane (Hij ~ 1 for i and j both in the transverse plane). This result means that 
fractality in multiparticle production is self-affine rather than self-similar. 

In order to show the independence of this conclusion of the particular set of variables 
being used, the same analysis has also been done with the set (y,pty,Ptz) instead of (y,pt, </?)• 
The corresponding results for the Hurst exponents are: 



/r " = TIT = °' 47 ± °- 06 



^p tz = ^ L = 0-47 ±0.06; 

*W = ^| = 0-99 ± 0.01 . 

They show that the rule H y j ps 0.5, « 1 for i and j denoting variables in the transverse 
plane also holds for the variable set (y,Pt y ,Ptz)- 



6 Self-affine analysis 

With the Hurst exponents obtained above, we can perform a self-affine analysis in three- 
and two-dimensional phase space. For convenience, we approximate the Hurst exponent for 
(y,p t ) and (y,ip) by 

Hyj =2> ti=Pt,<p), 

but use 

H ptl p = 1. 

From (3), it follows that A Pt = \ v = Xy. For p t and <p, we use a partitioning M y = 
1,2,3, ...,10 and M pt = M v = 1, 4, 9, . . . , 100 in two-dimensional analysis, but M y = 
1, 2, 3, . . . , 7 and M Pt = = 1, 4, 9, . . . , 49 in three-dimensional analysis. The results of 
the three-dimensional self-affine analysis on F2 are given by solid circles in Fig. 3. Those of 



the corresponding self-similar analysis are repeated by open circles, for comparison. In order 
to show the quality of the scaling law, linear fits 

\nF q = A + cj) q \nM y (7) 

are compared to the data in Fig. 3. The fit results are given in Table 2. To reduce the 
influence of momentum conservation, the first point is not used in the fits. 

Contributions to % 2 as shown in the top part of Table 2 are only from the diagonal terms 
of the covariance matrix. In fact, the diagonal terms provide the main contribution, and the 
relative size of x 2 for self-similar and self-affine analyses should not change dramatically by 
adding the contribution from non-diagonal terms that account for the correlation between 
points at different bin size. As is shown by comparing the results from fits to the unweighted 
sample, with and without considering the non-diagonal terms (two lower parts of Table 2), 
the fit results are retained better for the self-affine than for the self-similar analysis. A similar 
conclusion can be drawn from an inspection of Fig. 3, itself. While the self-similar analysis 
leads to an upward bending, this effect is absent in the self-affine analysis. Of course, higher 
statistics data would be needed to definitively prove this point. 

Even though the errors in the self-affine analysis are large, the results shown in Fig. 3 
support the expectation [7] that a 3-dimensional self-affine analysis would lead to the full 
increase of \nF q with increasing InM right from the beginning. On the other hand, in a 
self-similar analysis, the full increase would only be reached for large M, so that an up- 
ward bending would be observed. This upward bending is indeed present for the self-similar 
analysis. If future experiments can confirm the linear increase of the self-affine results with 
improved statistics, this would mean that the scaling law (1) is better observed in self-affine 
analysis than in self-similar analysis. 

The two-dimensional self-affine projections are presented in Fig. 4. It can be seen that 
lni*2 increases smoothly with increasing lnM y and the trends are similar for (y,pt) and (y, 92) 
(neglecting the first point in (y,ip)), meaning that in a self-affine analysis, the influence of 
artificial projection effects is reduced with respect to that observed in self-similar analysis. 

7 Conclusions 

In this paper we present a self-affine analysis of factorial moments in three-, as well as in 
two-dimensional phase space, on the NA22 data for 7r + p and K + p collisions at 250 GeV/c. 
The results are compared with those from a corresponding self-similar analysis. 

From fitting the factorial moments of the one-dimensional projections by the Ochs satu- 
ration formula (5), the Hurst exponents Hij are derived for all combinations of phase-space 
variables used. The value of for a combination of rapidity y with a transverse direction is 
approximately equal to 0.5. We conclude, therefore, that fractality in multiparticle produc- 
tion is self-affine, rather than self-similar. In the transverse plane, stays approximately 
equal to 1.0 and, therefore, shows merely self-similar fluctuation within that particular plane. 
Such a behavior can be understood from the privileged role of the longitudinal axis in mul- 
tihadron production and the symmetry within the plane transverse to this direction. This 
important point has been neglected in fluctuation analysis of multiparticle final states, so far. 

Furthermore, the three-dimensional self-affine analysis shows a better scaling behavior 
than does the corresponding self-similar analysis. The two-dimensional self-affine projections 
with = 0.5, i.e., (y,pt) and (y, if) turn out to show a behavior more similar to each other 
than in the corresponding self-similar projections. 



It would be interesting to see how the value of Hij changes with the type of collision and 
with incident energy. 
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Table 1. The parameter values obtained from a fit by (5). 



variables 


a 


b 


c 


Omitting 










point(s) 


y 


1.336 ± 0.005 


0.218 ± 0.042 


1.140 ±0.245 


1 


Pt 


1.534 ±0.021 


0.340 ± 0.021 


0.021 ± 0.006 


1 


f 


1.497 ±0.019 


0.420 ± 0.019 


0.014 ±0.005 


1 - 3 



Table 2. The parameter values obtained from a fit by (7). 





without bin-size correlation 


method 


A 


4>2 


X'/NDF 


weighted self-afhne 


-0.04 ± 0.03 


0.32 ±0.03 


7/4 


weighted self-similar 


0.11 ±0.01 


0.12 ±0.01 


9/4 




without bin-size correlation 


method 


A 




X 2 /NDF 


unweighted self-affine 


-0.08 ±0.02 


0.33 ± 0.03 


12/4 


unweighted self-similar 


0.09 ±0.01 


0.10 ±0.01 


20/4 




with bin-size correlation 


method 


A 


4>2 


X 2 /NDF 


unweighted self-affine 


-0.08 ±0.02 


0.34 ±0.02 


14/4 


unweighted self-similar 


0.10 ±0.01 


0.10 ±0.01 


32/4 



FIGURE CAPTIONS 



Fig. 1 Self-similar analysis of F q in the set of variables (y,pt, <p) in one, two, and three dimen- 
sions, as indicated. 

Fig. 2 Saturation curves for F2 in the three one-dimensional variables indicated. The curves 
are fits by (5) after omission of the first point (first 3 points in the case of -F2 (</?))• 

Fig. 3 Comparison of the three-dimensional self-affine and self-similar analyses. 

Fig. 4 Comparison of the two-dimensional projections of self-affine and self-similar analyses. 
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